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Long-range magnetic ordering in ifeBaNiOs (_R=magnetic 
rare earth) quasi- 1-dimensional mixed-spin antiferromagnets 
is described by a simple mean-field model that is based on the 
intrinsic staggered magnetization function of isolated Haldane 
spin chains for the Ni-subsystem, and single-ion magnetiza- 
tion functions for the rare earth ions. The model is applied 
to new experimental results obtained in powder diffraction 
experiments on Nd2BaNiOs and NdYBaNiOs, and to previ- 
ously published diffraction data for E^BaNiOs. From this 
analysis we extract the bare staggered magnetization curve 
for Haldane spin chains in these compounds. 
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The quantum-disordered ground state and the famous 
Haldane energy gap in the magnetic excitation spectrum 
[1] have kept one-dimensional (1-D) integer-spin Heisen- 
berg antiferromagnets (HAF) at the center of attention of 
condensed matter physicists for the last 15 years. Among 
the more recent developments are studies of such sys- 
tems in external uniform magnetic fields [2-4]. It was 
found that in sufficiently strong fields one of the three 
Haldane-gap modes undergoes a complete softening at 
some particular wave vector. The result is a transition 
to a new phase with long-range antiferromagnetic cor- 
relations (see for example Refs. [5,6]). The effect of a 
staggered field H ff , to which a Haldane chain is most 
susceptible, is expected to be no less dramatic. Unfortu- 
nately, this problem has been given much less attention 
in literature, simply because such conditions are almost 
impossible to realize in an experiment. The only chance 
of producing a magnetic field modulated on the micro- 
scopic scale is to make use of some periodic modulation 
that is intrinsic to the system under investigation. This 
can be, for example, a structural modulation, as is the 
case in NENP, one of the best-known Haldane-gap com- 
pounds (Refs. [7,8] and references therein). The S = 1 
chains in this material consist of alternating crystallo- 
graphically non-equivalent Ni 2+ ions with slightly differ- 
ent gyromagnetic ratios. A weak effective staggered field 
can thus be induced in NENP by applying a uniform ex- 
ternal field [9] . Unfortunately, the effect of the staggered 
component is obscured by the response of the system to 
the strong uniform field itself [9,10]. 

A more direct approach is to use an intrinsic magnetic 



modulation in a material that, in addition to integer-spin 
Heisenberg chains, has other magnetic ions. Should the 
latter become ordered magnetically with an appropriate 
propagation vector, they will project an effective stag- 
gered exchange field on the Haldane spin chains. The 
magnitude of the staggered field can be varied in an ex- 
periment indirectly, by changing the temperature, and 
thus the magnitude of the magnetic order parameter. By 
measuring the induced moment on the Haldane chains 
one can hope to directly measure the staggered magneti- 
zation curve for a Haldane spin system. This function, 
that we shall denote as A / [(i? 7r ), is one of the princi- 
pal characteristics of a quantum-disordered integer-spin 
HAF, but to date has not been determined experimen- 
tally. In the present paper we shall demonstrate how 
Al(-ff^) can be determined by properly analyzing the re- 
sults of simple neutron diffraction experiments. 

It appears that the best experimental realizations of 
the mechanism described in the previous paragraph are 
to be found among linear-chain nickelates with the gen- 
eral formula i?2BaNiC>5 (i?=magnetic rare earth). The 
5 = 1 Ni 2+ ions in these compounds are arranged in dis- 
tinct chains that run along the a axis of the orthorhom- 
bic crystal structure. In-chain interactions between the 
Ni spins are antiferromagnetic, with the exchange con- 
stant J rs 300 K, and interchain coupling is negligi- 
ble. The magnetic i? 3+ ions are positioned in between 
the Ni chains, and are expected to be weakly coupled 
to the Ni 2+ moments (exchange parameters of the or- 
der of tens of Kelvin) and even weaker between them- 
selves (coupling constants of less than 1 K) [11]. All 
i?2BaNiC>5 species order antiferromagnetically with Neel 
temperatures ranging from 16 to 80 K [11,12]. In the 
ordered phase a non-zero ordered moment is found on 
both R 3+ and Ni 2+ sites. The magnetic order parame- 
ter of Ni 2+ never achieves complete saturation, and at 
T — ► is substantially smaller than the expected value 
of 2/ib per ion, a hint at that quantum spin fluctua- 
tions in the Ni-chains survive well below the Neel tem- 
perature Tn. Haldane gap excitations propagating on 
the S = 1 Ni-chains have been observed both above 
and below the Neel temperature in all materials studied 
so far, including Pr 2 BaNi0 5 [13], Nd 2 BaNi0 5 [14], and 
(Nd a; Yi_ 2 ;)2BaNi05 [15]. These spin excitations have a 
purely 1-D dispersion and dynamic structure factor, and 
are practically indistinguishable from Haldane modes in 
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Y 2 BaNi05, an extensively studied reference Haldane gap 
system, where no long-range magnetic ordering occurs 
even at low temperatuures [16-19]. The basis of our 
current understanding of i? 2 BaNi05 compounds is that 
the Haldane singlet ground state of individual S = 1 
Ni-chains is preserved in the magnetically ordered phase 
[13,20]. The non-zero ordered moment on the Ni 2+ sites 
is viewed as a result of polarization of the quantum- 
disordered 5 = 1 chains by an effective staggered ex- 
change field, generated by the ordered i?-sublattice. 

Based on this physical picture, we shall construct a 
simple mean-field (MF) model of magnetic ordering in 
i? 2 BaNiOs materials. In the conventional MF approach 
all magnetic ions are treated as isolated moments that 
become polarized by the effective exchange field. In the 
ii^BaNiCVs however, in-chain Ni-Ni exchange coupling 
is dominant over all other magnetic interactions. It is 
therefore more appropriate to treat only the rare earths 
ions as individual moments [11]. In the spirit of what 
is said above, the bare magnetization curve for the Ni- 
sublattice should be taken in the form of the staggered 
magnetization function M. for an isolated S = 1 quan- 
tum spin chain. Ni-Ni interactions are not to be included 
explicitly into the MF theory, since they are already built 
into the function M.. The MF equations for the Ni- 
subsystem can thus be written in the following form: 

M< Ni ) = gSnBM(H^), (1) 
=2a R M( R \ (2) 

Here M < Ni ) and are the sublattice magnetizations 

(per ion) for the Ni and R subsystems, respectively, i^( Nl ) 
is the effective staggered field that acts on the Ni-chains, 
g = 2 is the gyromagnetic ratio for Ni 2+ , S = 1 is the spin 
of Ni 2+ , and a R is the effective MF coupling constant. 
To complete the model, we have to write down the MF 
equations for the i?-subsystem as well. Unfortunately, 
due to the low site-symmetry for R 3+ ions in i? 2 BaNi05 
materials [13], modeling the exact bare magnetization 
curve of the rare earths is impossible without knowing 
the details of their electronic structure. Nevertheless, 
for those compounds in which R 3+ is a Kramers ion, we 
can hope to get a reasonable approximation using the 
Brillouin function: 

M« = M^ R) tanh(# W M (R) /«T), (3) 
= a R M<- Ni >. (4) 

( B) 

Here Mq is the effective moment of the rare earth ion 
and H( R ~) is the mean field acting on the rare earths. 
Direct R-R interactions are expected to be weak, and 
Tjy is predominantly defined by the magnitude of i?-Ni 
exchange. This fact enables us to simplify the equations 
by omitting R-R coupling in Eq. 4. 

Our "semi-quantum" model differs significantly from 
the usual MF model for classical magnets. In the latter, 



at T — > all sublattices become fully saturated as the 
bare single-ion susceptibility diverges as 1/T. In con- 
trast, in our model the bare staggered susceptibility of 
the Haldane chains remains finite at T = 0. Moreover, 
at T < A (A « 9 mcV, or « 105 K i? 2 BaNi0 5 com- 
pounds), the function M.{H^) is expected to be almost 
T- independent. Our model will therefore produce quali- 
tatively different results for the T-dependencies and sat- 
uration values of the magnetic order parameter for Ni. 

It has to be noted that the idea of treating interchain 
coupling at the MF or RPA (Random Phase Approxi- 
mation) level is, in itself, not a new idea, but a well- 
established technique [21]. In their pioneering work on 
CsNiCIa, the first Haldane-gap material studied experi- 
mentally, Buyers et al. [22] , and Affleck [23] implemented 
this approach to explain magnetic ordering and calculate 
the spin wave dispersion relations. The main difference 
between CsNiCl3 and i? 2 BaNi05 is that exchange cou- 
pling between individual Haldane chains is direct in the 
former system, and mediated by the rare earth ions in 
the latter. For directly coupled Haldane spin chains, the 
magnitude of interchain interactions must exceed some 
critical value in order for the system to order magneti- 
cally [23]. In a MF treatment, the temperature depen- 
dence of the ordered moment is defined by the intrinsic 
temperature dependence of the susceptibility of individ- 
ual chains. In contrast, in our case of i? 2 BaNi05 com- 
pounds, magnetic ordering is driven by the 1/T-divergent 
susceptibility of the rare earth subsystem. At sufficiently 
low temperature long-range order will therefore occur 
for arbitrary small i?-Ni interactions, and in the case of 
TV A, we can use the approximation where the bare 
magnetization curve of isolated chains is T-independent. 
i? 2 BaNi05 materials are thus more convenient as model 
systems for studying the effect of staggered field on Hal- 
dane spin chains. 

To test the validity of our model we shall compare its 
predictions for the temperature dependence of sublattice 
magnetizations with experimental results obtained using 
neutron diffraction. Since one of our primary goals is to 
verify that the function M. is indeed T-indepcndent at 
T < A, we need to compare experimental M(T) curves 
measured in material with substantially different order- 
ing temperatures. Moreover, as the constant a R is clearly 
dependent of R, the only way to exclude T as an im- 
plicit variable is to study several dilute systems with the 
formula (i? a; Yi_ a ;) 2 BaNi05. In such compounds part 
of the magnetic R 3+, s are randomly replaced by non- 
magnetic Y 3+ ions. Tjv is strongly dependent on the 
rare earth concentration x [15]. At the same time, we 
can expect the coupling constant a R to be the same in 
all (ii x Y 1 _ a; )2BaNi05 systems, provided in Eq. 2 
is replaced by the average moment on the R sites, i.e., 
by = xM( R \ Indeed, in the crystal structure each 

Ni site is coordinated to four i?-sites, and any effect of 
disorder, for not too small i?-content, is reduced by aver- 



2 



aging over nearest neighbors, and the effective staggered 
field scales proportionately to x. 

To date, accurate T-dependent data exist only for 
Ho 2 BaNi0 5 (T N = 53 K) [24] and Er 2 BaNi0 5 (T N = 
33 K) [25]. Of these materials only the Er nickelate fits 
the condition of having a Kramers rare earth ion, essen- 
tial for using Eq. 3. To have more compounds to work 
with, and to apply our model to systems with substan- 
tially different ^-concentrations x, we have performed 
new powder neutron diffraction studies of Nd 2 BaNi05 
(T N = 48 K, x = 1) and NdYBaNiOg {T N = 29 K, 
x = 0.5). The experiments were done on the posi- 
tion sensitive detector diffractometer DIB at the Insti- 
tut Laue-Langevein, Grenoble, France. Details of the 
rather straightforward experiment and data analysis will 
be reported elsewhere. At all temperatures the obtained 
magnetic structures in both materials are defined by the 
propagation vector (1/2,0,1/2), and are similar to that 
in Ho 2 BaNi0 5 [24]. In NdYBaNi0 5 the Ni 2 + and Nd 3 + 
moments lie in the (a, c) crystallographic plane, and form 
almost temperature-independent angles of </>( Nl ) « 35° 
and </>( Nd ) w 0° with the c-axis, respectively. The same 
applies to Nd 2 BaNi0 5 at T < 40 K. Between 40 and 
48 K however, i. e., just below T/v, both the Ni and Nd 
ordered moments in Nd 2 BaNiC>5 undergo a significant 
re-orientation, in agreement with Rcf. [12]. Our princi- 
pal experimental results, the temperature dependencies 
of the ordered moments, are summarized in Fig. 1. In 
the same figure we show the data for Er 2 BaNiOs, taken 
from Ref. [25]. In Er 2 BaNiOs the orientation of the or- 
dered moments of Ni 2+ and Er 3+ is T- independent, with 
(f) (N V = 65° and ^ = 88°, respectively. In the Ho- 
system the orientation of ordered moments change only 
slightly with decreasing temperature, and the angles are, 
correspondingly, 0< Ni ) w 25° and 0< Ho ) w 0° [24]. 

We are now in a position to plot M^ Nl ) as a function 
of M (i?) . For Nd 2 BaNi0 5 and NdYBaNiOg (x = 1 and 
x = 0.5, respectively) this curve is shown in Fig. 2 in 
solid triangles and squares, correspondingly. For these 
two materials in which the ordering temperatures differ 
by almost a factor of 2 we obtain an excellent data col- 
lapse, which confirms that the function M. is almost T- 
independent. To place the Er 2 BaNi05 data on the same 
plot we have to re-scale the absciss to compensate for 
the difference between a^d and «Er- This is done in such 
a way as to obtain the best overlap between the three 
curves. The best data collapse is achieved by scaling 
down the Er moment by a factor of 3.2(1) ( Fig. 2, open 
circles). For Ho 2 BaNi05, a non-Kramers system, the ap- 
proximation (3) is not appropriate, yet Eqs. (1,2) are ex- 
pected to remain valid. We can therefore plot the existing 
data for Ho 2 BaNi05 [24] in the graph on Fig. 2, applying 
the same absciss- rescaling procedure as for Er 2 BaNi05. 
The best data collapse in this case is obtained by scaling 
down the Ho moment by a factor 4.2(1) ( Fig. 2, open di- 
amonds). We see that, despite the different -R-substitute, 



the overall shape of the M^ Nl ^ vs. curve in all four 

materials is very similar, in support of our conjecture 
that M is an intrinsic property of the Ni chains and is 
the same for all i? 2 BaNi05 systems. For further use we 
shall fit the cumulative data in Fig. 2 with a purely em- 
pirical function: 

gSn B M{a R M^) = Aarctan^M^), (5) 

A good fit (solid line in Fig. 2) is obtained with A = 
1.17(4) n B , £>Nd = 1.7(1) \i B \ B Er = £? Nd /3.2 = 
0.43(3) fi B \ and B Ho = B Nd /4.2 = 0.33(3) fi B \ 

In the context of our model the curve in Fig. 2 is noth- 
ing else but the staggered magnetization function for a 
Haldane spin chain that we are trying to determine in 
this study. The only thing we are missing at this point 
is a proper scale (units of magnetic field) on the absciss, 
which is currently labeled in units of magnetization. Con- 
version to field units is done by the coupling constant 
a, that we shall obtain by analyzing the temperature- 
dependent experimental data with our mean-field equa- 
tions. The analysis can be performed for R =Nd/Y and 
R =Er, but not for R =Ho, where, as discussed above, 
Eq. 3 does not apply. In place of Mq R ^ in Eq. 4 we 
shall use the saturation moment for the rare earth ions 
in each system. That the values of Mq N<1 ^ are slightly 
different in the two Nd-based compounds (2.68/ib and 
2. 3^b, respectively) only emphasizes the fact that Eq. 4 is 
no more than an approximation to the actual magnetiza- 
tion function for an isolated R 3+ ion. Combining our MF 
equations with Eq. 5 we use least-squares refinement to 
fit the experimental data in Fig. 1, having a_R as the only 
adjustable parameter for each compound. Very good fits 
are obtained for all systems (solid lines in Fig. 1). The 
only discrepancy is above 40 K for Nd 2 BaNi05, where the 
magnetic structure changes dramatically and our model 
is not applicable anyway. The refined values for the 
coupling constants are: a^j = 0.82(2) • 10 5 Oc/^s in 
Nd 2 BaNi0 5 , a$} = 0.99(2) • 10 5 Oc/^ B in NdYBaNiOg, 
and a Er = 0.183(3) • 10 5 Oc/fi B in Er 2 BaNi0 5 . Obtain- 
ing similar values of aNd f° r the tw° Nd-based nicke- 
lates is another self-consistency check for our approach. 
Equation 2 can now be used to replot the induced mo- 
ment on the Ni-chains as a function of effective exchange 
field. As the data collapse in Fig. 2 is rather good, 
for this purpose we can use the effective average value 
«Nd = {a$ A + a$ A + 3.2a E r)/3 = 0.8(1) • 10 5 Oe/^ B . 
The resulting staggered field scale is shown in the top 
axis in Fig. 2. 

Our model, simple-minded as it is, appears to ade- 
quately describe the magnetic ordering in all i? 2 BaNi05 
systems studied so far. Its application provides the first 
measurement of the staggered magnetization curve for a 
Haldane spin chain. We have recently succeeded in calcu- 
lating this function theoretically, extending the approach 
laid out in Ref. [20]. This calculation is in good agree- 
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ment with our experimental data, but the discussion is 
beyond the scope of this paper and will be addressed in 
a forthcoming publication. 
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FIG. 1. Measured temperature dependencies of the Ni- 
(open symbols) and R- (solid symbols) ordered moments in 
Nd 2 BaNi0 5 and NdYBaNiOs (this work), and Er 2 BaNiO s 
(digitized from Ref. [25]). The solid lines are single-parameter 
fits using the mean field model described in the text. 

FIG. 2. Ordered staggered moment on the Ni-chains plot- 
ted against the ordered moment of the _R-sublattice (bottom 
axes). The solid line is an empirical fit, as described in the 
text. The data for Ho2BaNiOs and E^BaNiOs are taken from 
Refs. [24,25]. The top axis shows the corresponding staggered 
exchange field acting on the Ni-chains. 
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